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We extend Multichannel Quantum Defect Theory (MQDT) to ultracold collisions involving high 
partial wave quantum numbers L. This requires a careful standardization of the MQDT reference 
wave functions at long range to ensure their linear independence. To illustrate the simplicity and 
accuracy of the resulting theory, we perform a comprehensive calculation of L < 2 Fano-Feshbach 
resonances in the range — 1000 G for the scattering of ^"K + *^Rb in their lowest hyperfine states. 



I. INTRODUCTION 

The constituents of cold gaseous matter continue to 
grow in variety and complexity far beyond their origins 
in alkali-metal atoms to encompass open shell atoms, 
molecules, free radicals, and ions [1]. A detailed under- 
standing of their collision processes is crucial in deter- 
mining the properties and prospects for control of these 
gases. 

However, along with the growing complexity of the 
atoms and molecules involved, the difficulty of accu- 
rate scattering calculations grows as well. High-spin and 
open-shell atoms contain a multiplicity of internal states, 
and molecules incorporate rotational and vibrational de- 
grees of freedom. At ultralow temperatures all these 
degrees of freedom must be accounted for because they 
all describe energies that are typically large compared 
to translational kinetic energies in the gas. Moreover, 
anisotropic species are likely to involve angular momen- 
tum partial waves far larger than the single value L = 
that often dominates alkali atom cold collisions. These 
considerations can lead to enormous complexity, even in 
seemingly straightforward problems. Consider for ex- 
ample the scattering of Li atoms with ground state Li2 
molecules, which necessitates the simultaneous solution 
of thousands of coupled Schrodinger equations 

It is therefore worthwhile to explore alternative meth- 
ods of scattering theory that are less computationally in- 
tensive, yet still accurate. A particularly appealing can- 
didate for this purpose is the Multichannel Quantum De- 
fect Theory (MQDT), which was originally developed to 
describe complex spectra of atoms [3] [1] [S] . It has more 
recently seen fruitful application to atomic collisions ^ 
and, in particular, cold collisions [7] [H] [9]. 

In recent years various approaches have been applied 
to ultracold collisions with an aim of describing and pre- 
dicting scattering observables and, in particular, the lo- 
cations and widths of Fano-Feshbach resonances. For ex- 
ample, recent developments in MQDT and related ideas 
have stressed the simplicity and analytic behavior of the 
theory [lQ][n][l2][l3][n![l5][l6][l7][l8]. Also, the 
analysis of long-range s-wave solutions developed in [T5] 
adopts semiclassical ideas and creates a good approxi- 
mation to the long-range field solutions that resembles 
quantum defect ideas to a degree. To a good approx- 
imation, three parameters can describe cold collisions 



of alkali-metal atoms. There are two scattering lengths 
that describe the physics at small interparticle separa- 
tion R and one dispersion coefficient Cg that describes 
the physics at large R [H]. In parallel developments, the 
asymptotic bound state method stresses direct numeri- 
cal diagonalization in a basis of singlet and triplet states 
that are coupled by hyperfine and magnetic interactions 
PP] [2T] [55]. This has generated extremely rapid and 
accurate numerical calculations in these cases. 

The present paper describes an attempt to retain the 
analytic structure of MQDT and simultaneously perform 
numerically exact scattering calculations. On the one 
hand, MQDT allows the "hard" short-range physics to 
be described via parameters that are nearly energy- and 
field-independent. In addition, analytic expressions are 
derived for some of the long-range parameters in the 
threshold limit. On the other hand, our treatment is 
numerically exact, even in the case of realistic long-range 
potentials that are not purely of — Cg/i?^ character and 
that possess high centrifugal angular momentum. As 
a numerical method, it is also fast because it eschews 
repeated numerical calculation of the short-range wave 
fimctions. 

MQDT takes into account the natural separation of 
length and energy scales inherent in collision problems. 
Namely, the sensitive dependence of collision observables 
on energy and electric or magnetic fields arises from in- 
teractions between the scattering partners at large R. 
These dependencies are described ultimately through a 
carefully chosen set of reference wave functions in this re- 
gion. This part of the calculation is relatively fast, as the 
reference wave functions are defined separately for each 
set of quantum numbers identifying a scattering chan- 
nel. Vice versa, the "hard" part of the calculation lies at 
smaller R, where the channels are strongly coupled to- 
gether and must all be dealt with at once. However, in 
this circumstance, the energy scales driving the physics 
are far larger than the - mK scale of cold collisions. A 
properly chosen representation of the small- i? wave func- 
tion can then be quite weakly energy dependent, allow- 
ing ready interpolation in energy and field that greatly 
reduces computational time 10 . 

However, accuracy of the MQDT method relies on the 
availability of a pair of accurate, linearly independent ref- 
erence wave functions in each channel. This can become 
problematic in cases of high partial wave angular mo- 
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mentum and extremely low collision energy, where there 
exists a substantial region of classically forbidden motion 
underneath a centrifugal barrier. Such a barrier is prob- 
lematic because two reference wave functions that are 
perfectly linearly independent prior to entering the bar- 
rier can often become approximately linearly dependent 
(hence useless to the theory) underneath and beyond this 
barrier. 

Our main focus here is to determine reference wave 
functions in a way that maximizes their linear indepen- 
dence, even underneath the centrifugal barrier. The 
method is then applied to the scattering of potassium 
and rubidium atoms, for which a reasonable and ex- 
perimentally constructed Hamiltonian already exists and 
the long-range interaction includes contributions from Cg 
and Cio in addition to the usual Cg. MQDT reproduces 
the full close-coupling calculation quite accurately with 
orders of magnitude less computational effort. In par- 
ticular, this version of MQDT reproduces magnetic field 
Fano-Feshbach resonances accurately, even those that re- 
side in high partial wave states. 



II. THEORY 

A. Scattering in the MQDT Picture 

In two-body atomic or molecular scattering, the wave 
function can be expanded in a basis of N magnetic-field- 
dressed hyperfine channels that include the partial wave 
quantum number L^, 



JV 



V' = i?-i^<i>.(n)V',(i?) 
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1=1 



where f2 represents all angular coordinates and spin de- 
grees of freedom. This wave function satisfies a set of 
coupled radial Schrodinger equations involving the po- 
tential matrix V{R), 
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Ei is the channel collision energy, i.e., the energy above 
threshold for each channel Ei — E — £:^hrosh^ where E is 
the collision energy and £;thrcsh ^-j^g threshold energy 
of the i-th channel. 

Here and throughout this paper (unless otherwise spec- 
ified), all lengths are in units of the natural length scale /3 
of the potential V , and all energies are in units of the nat- 
ural energy scale Ep = It? /2^[3'^, where ^ is the reduced 
mass. As a consequence, ^ is scaled out of many equa- 
tions, and ki = \/Ei — Vij has units of I3~^. The form of 
V at long range suggests the value of /3. For example, if V 
is well described at long range by the isotropic potential 
-Ce/R^, the natural unit of length is /3 = {2nCe/h'^)^/'^. 
Later sections discuss the collision of ""^K + ®^Rb, for 



which /3 = 143.9 slq and Ep = 152.7 /iK, where ao is the 
Bohr radius. 

One can calculate scattering observables by solving 
equation ^ subject to physical boundary conditions. In 
each channel tp the boundary conditions ipi — a,t 
R ^ 0. The closed-channel {Ei < 0) components of ip 
vanish at R = oo. Hence, if Ng is the number of open 
channels {Ei > 0), there exist No independent ip's which 
only have a non- vanishing amplitude in the open channels 
asymptotically. The open-channel components of these 
wave functions are energy-normalized and represent the 
No X No asymptotic solution matrix 4". 

Since the potential goes to zero in the limit i? — >■ oo, 
the open-channel wave functions become linear combina- 
tions of sine and cosine, and ^' can be written in terms 
of energy-normalized, free-particle solutions. 



i^ihR) (jLAhR)^, - nL,{hR)K- 
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where and ul are spherical Bessel functions of the first 
and second kind, respectively, 

X) sm{kiR — LiTT /2) 

~^ hR ' 
X) cos{kiR - LiTT/2) 



(4a) 
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Here, 



Sij is the Kronecker delta function, and ki — 

' V^i > 



Ei. Equation i3h defines the physical 



iiT-matrix which contains all the information necessary 
to compute scattering observables, including resonance 
behavior and threshold effects, and is simply related to 
the scattering matrix, 

S-Phys ^^ j^ .^phys -J _ ^^phys ^ - 1 ^ 

where / is the identity matrix. 

For many problems, such as atoms and molecules in- 
teracting via van der Waals potentials, the channels are 
approximately uncoupled beyond a radius R^^- In gen- 
eral, there are N independent solutions to equation ^ 
that have only the boundary conditions tpi = at R — 
and represent the N x N solution matrix M. Hence, 
matching M to single-channel reference wave functions / 
and g at R = R^ defines a short-range iiT-matrix [23| , 



(6) 



Here, / and g are solutions to the uncoupled radial 
Schrodinger equations in the long-range potential V^'^, 

U{U 

'di?2 + 



1) 



i?2 



V " - E, 



9i 



(7) 



The matching is best done when all channels are locally 
open {Ei > Vij{Rni)) because / and g are oscillatory and 
can, therefore, easily be made linearly independent. To 
the extent that the channel coupling is negligible beyond 
-Rnu a-pplying boundary conditions at i? = oo in terms 
of / and g allows scattering observables to be computed 
accurately. 
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B. Reference Wave Functions 

The major benefit of using tlie reference wave func- 
tions / and g is that they do not need to satisfy physical 
boundary conditions. Their boundary conditions can in- 
stead be smooth, analytic functions of energy. In partic- 
ular, choosing WKB boundary conditions well within the 
classically allowed region at i? = i?x < accomplishes 
this goal d] [H], 



1 r'^ 

/'(^) = /, sin( / h{R')AR' + at i? = i?, , 

yJki{R) Jr^ 



= -^r7^^°^( rhiR')dR' 



at i? — -Rx 7 
(8b) 



where can be any channel-dependent phase that is con- 
stant in R and energy [10]. The set of equations Q and 
their full radial derivatives define / and g. As the WKB 
boundary conditions define these single-channel reference 
wave functions at the single radius i?xi they still consti- 
tute exact solutions of the radial Schrodinger equation 

Moreover, the particular solutions defined by (|8| have 
several advantages. First, unlike the wave functions Q, 
these reference wave functions are not energy-normalized. 
They are well defined for Ei < and analytic in en- 
ergy across threshold. Second, as the large kinetic energy 
at short range dominates the low collision energies and 
relatively small Zceman shifts of typical cold collisions, 
WKB boundary conditions lead to reference wave func- 
tions that are weakly dependent on collision energy and 
magnetic field. From this follows the weak energy and 
field dependence of K^^' . Third, this particular choice 
of boundary conditions allows / and g to be maximally 
linearly independent at short range. 

The matrix K^"^ and the linearly independent refer- 
ence wave functions / and g contain all the information 
necessary to calculate scattering observables. The quan- 
tum defect theory of [HI defines the four parameters rj, 
A, Q, and /Jsurko that describe the asymptotic behavior 
of the wave functions / and g. Hence, these parame- 
ters are also smooth functions of collision energy and 
magnetic field and completely describe the long-range 
physics. The notation in this paper only differs from the 
notation of JIO] by the use of 7 instead of /3Buikc, where 
cot 7 = tan/3Burko- The introduction of 7 emphasizes 
its relationship with Q , and section [TV] demonstrates this 
relationship. 

The calculation of 5^^^** requires two linearly indepen- 
dent, energy-normalized wave functions at large R in each 
open channel and the bound state wave function in each 
closed channel. To this end, the parameters A and Q 
create a Wronskian-preserving transformation between 
the reference wave functions / and g and two energy- 



normalized wave functions / and g for Ei > Q, 
f\_( A-l^ ^//^ 



(9) 



The parameter A is responsible for the energy- 
normalization of / and (7, and the parameter Q accounts 
for the different phase g accumulates in y'"' than does /. 
The phase shift 77 describes how / and g differ from the 
spherical Bessel functions asymptotically, 

/j "•^°°> ^ sin(fc.,i? - Li-Kjl + m) , (10a) 
gi ^■^°°) --^ cos{hR- Lin/2 + r]i) . (10b) 

For < 0, the parameter 7 determines the hnear com- 
bmation of / and g that vanishes as i? — ^ 00, 



tan-fjt+gi ^cx e '^'^ 



(11) 



where Ki — iki. 

Calculating the four MQDT parameters requires the 
evaluation of several Wronskians that involve / and g as 
R — >■ 00, 
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(12c) 



(12d) 



R^ca 



where W (2/1,2/2) is the Wronskian with respect to R of 
any two functions j/i and 2/2 1 

W{yi,y2) = yi{R)—^ 2/2 (i?) ■ (13) 

The MQDT parameters directly translate K^'^ into ob- 
servables. By partitioning K^^ into open (P) and closed 
(Q) channels, simple algebra produces the physical scat- 
tering matrix [lOj . 



K = Kfp - K^'^{K^4^ + cot7)-^i^§p 

K = A^/^K{i + gky^A^'^ , 



^phys _ gi)) 



I + iK 
I-iK^ 



(14a) 
(14b) 

(14c) 



4 




mfS/ 

-1 

-2 

0.4 0.8 1.2 

(a)/ asymptotically coincides with x+ at zero energy. 
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(b)g asymptotically coincides with x- at zero energy 



FIG. 1: (Color online) The reference wave functions / (red 
dashed curve) and g (green dotted curve) are shown for 
Ei « 6.547 X 10"^ in the long-range potential F'" = ~1/R^ 
(blue solid curve) with L = 2. At short range, / and g are lin- 
early independent. Under the classical barrier, these functions 
(a) lose their independence when / is chosen to asymptotically 
coincide with x+ a-t zero energy and (b) retain their indepen- 
dence when g is chosen to asymptotically coincide with X- at 
zero energy. 



III. STANDARDIZING MQDT 

Since the MQDT parameters clearly depend on the 
particular choice of reference wave functions, standardiz- 
ing this choice allows the MQDT parameters for a partic- 
ular long-range potential to be tabulated once and for all 
and defines a simple procedure to find K^^ [TU]. In gen- 
eral, the boundary conditions (|8| define an infinite family 
of reference wave functions - one set for each value of 0^. 
This section identifies a value of (j)i that guarantees both 
the maximal numerical stability of the MQDT parame- 
ters as well as their smooth, analytic energy behavior. 
Moreover, since the calculation of (j)i is required to ap- 
ply the boundary conditions (8]), this calculation must be 
numerically stable - even at high L - to be useful. 



In the manner of [TU] , the MQDT parameters are stan- 
dardized by letting one of the reference wave functions 
asymptotically coincide (up to a normalization) with a 
particular wave function at zero energy. The choice of 
this standard, zero-energy wave function, therefore, iden- 
tifies a particular value of (pi and determines a particular 
/ and g at zero energy. The energy dependence of / and 
g near threshold is then smoothly obtained from their 
WKB boundary conditions ([8|. In principle, any value 
of (pi is equally valid as long as / and g remain linearly 
independent. However, at ultracold energies, motion un- 
der a classical centrifugal barrier can often cause a pair 
of reference wave functions to become linearly dependent 
numerically, so it is our goal to find the value of that 
leads (in the limit of zero energy) to maximally linearly 
independent reference wave functions under and beyond 
this barrier. 

For all potentials that fall off faster than 
the asymptotic linear combination of and 
uniquely identifies a particular zero-energy wave func- 
tion. Consequently, the asymptotic linear combination 
of the analytic zero-energy solutions [25 , 



X+ = A/i?J_i(2L+i)(l/2i?2) il^cx R 

= /RJi(2i+i)(l/2i?2) i^cx R 



L+l 
-L 



(15a) 
(15b) 



identifies a particular zero-energy wave function in any 
2 potential that is dominated by V^'^ = —1/R^ asymptoti- 
cally. Here, J is the Bessel function of the first kind. For 
L > 0, these solutions take their asymptotic form under 
the classical barrier at R> 1. Therefore, all zero-energy 
wave functions start to resemble x+ in ttiis region - losing 
their linear independence numerically - except for those 
directly prop ortional to x_. 



Fig. l(a 



shows / and g in the long-range potential 
V'' = -LJW with Ei « 6.547 x 10"^^ ^his energy cor- 
responds to 1 fiK for -I- «^Rb when Cg = 4.300 x 10^ 
in atomic units. Here, / is chosen to asymptotically co- 
incide with x+ at zero energy. The boundary conditions 
^ ensure the maximal independence of / and g at short 
range, but they both resemble x+ in classically for- 
bidden region R > I. While / is required to grow like 
i?^"*"^ under the classical barrier, g also quickly begins to 
grow in a similar way. Hence, this choice for / leads to 
a set of reference wave functions that exhibit increasing 
linear dependence as the collision energy approaches zero 
and the classical barrier grows. In fact, g always asymp- 
totically diverges at zero energy except for a unique value 

of 04 . 

Since only x- remains numerically linearly indepen- 
dent from x+ in the limit i? — > c», letting g asymptoti- 
cally coincide with X- at zero energy guarantees that the 
zero-energy limit of g is maximally independent from / 
not only at short range but also well into the classically 
forbidden region. Moving away from zero energy causes 
g to gain a contribution from x+ asymptotically, but the 
classically forbidden region becomes smaller. Fig. 1(b) 
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shows that this choice for g leads to reference wave func- 
tions that are hnearly independent at both short range 
and long range - even at ultracold energies and high L. 
Hence, these reference wave functions are ideal for a nu- 
merical calculation of the MQDT parameters at ultracold 
energies. 

In order to implement this standardization, one must 
determine the value of (pi that is used to define / and 
g. For long-range potentials in which the zero-energy 
solutions are known analytically at all i?, this value of 
4>i is easily derived using the values of x+ and/or X- at 
_Rx- However, since even the zero-energy solutions are 
only known analytically for a limited number of power 
law potentials, calculating (pi numerically allows the use 
of the true long-range potential for a given scattering 
problem. 

To this end the zero-energy reference wave functions 
for (fii = Q can be numerically propagated from their 
boundary conditions at i?x << 1 to large R » 1. 
At i? >> 1 the reference wave functions are well ap- 
proximated by linear combinations of solutions that are 
known analytically, and tan 0^ is a simple ratio of two 
Wronskians. For example, the value of tan (pi that lets / 
asymptotically coincide with x+ at zero energy is given 

by, 



tan< 



W (x+,/(0,=o)) 

w {x+,g{4>,=Q)) 



(16a) 



An alternative choice allows g to asymptotically coincide 
with X- at zero energy, and this value of (j)i is given by. 



tanffli > 



W (x-,g(0.^Q)) 
W (x-Ji4>,=o)) 



(16b) 



The stable calculation of (pi requires the linear indepen- 
dence of the reference wave functions /(^^^o) and g(0.=o)- 
As they each represent a particular linear combination 
of x+ and X- asymptotically, their linear independence 
relies on their different, and quickly vanishing, contribu- 
tions from X-- Hence, they can easily lose their linear 
independence as R grows, causing the calculation of (pi 
to become unstable. For example. Fig. [2] demonstrates 
the difficulty of numerically calculating pi by showing 



the numerically computed ratio of Wronskians in ( 16a 



for all values of R. These calculations use the long-range 
potential V^' ^ -Cq/R^ - Ca/R^ - Cw/R^° for sev- 
eral values of L, where the dispersion coefficients are 
Cg =4.300x10^, Cg = 4.823xl0^ and Cio = 6.181x10^ 
in atomic units. These values are realistic for collisions 
of K Rb [ig. 

Fig. [2] shows that the calculation of tan pi , at least 
for i = 0, converges quickly after i? = 1, where the 
—Cs/R^ and — Cio/i?^° terms of the long-range poten- 
tial become dominated by the —Cq/R^ term. However, 



as 



L increases, a larger value of R is required to con- 
verge this calculation, and Fig. [2] shows that a converged 




FIG. 2: The convergence of tan^i with R is shown when 
L = — 3. Here, / is chosen to asymptotically coincide with 
X+ at zero energy in the long-range potential V^^ = —Ce/R^ — 
Cs/R^ - Cio/R^°, and 7£x = .1. 
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FIG. 3: The convergence of tan^i with R is shown when 
L — 0, 2, 10, and 20. Here, g is chosen to asymptotically 
coincide with X- at zero energy in the long-range potential 
V^' = -Ce/R'' - Cs/R^ - Cio/R^°, and R^ = .1. 



calculation of tan pi is not numerically stable for L > 1. 
Appendix|A]explores this instability and determines that, 
for L > 1, the asymptotic contribution to each reference 
wave function from x- is eventually dominated by the 
presence of finite numerical noise, causing the calcula- 
tion of tan^j to be unstable. If even (f>i cannot be stably 
computed, calculating the MQDT parameters is hope- 
less. 
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(12 1. Of particular importance is the parameter Q. The 





tan d 



FIG. 4: (Color online) The red dots are the MQDT pa- 
rameter Q shown as a function of tan<^i. The vertical lines 
represent the value of tan 4)i when (blue dashed line) / asymp- 
totically coincides with x+ ^t zero energy and (green solid 
line) g asymptotically coincides with x- ^ero energy. 
For this calculation, V'"' = -Ce/R^ - Cs/R^ - Cio/R'", 
where the dispersion coefficients are appropriate for K -f- Rb; 
Ei « 6.547 X 10"'', which corresponds to 100 nK; L = 0; and 

iix = .1. 



However, Appendix |A] also shows that this instability 
vanishes if g is chosen to asymptotically coincide with 
X- at zero energy - the same choice as in Fig. |l(b) For 
only this choice of g, the calculation of does not require 
finding the asymptotic contribution to /(0.=o) and g^^.^Q) 
from x_. Hence, in this case, the calculation of tancjii is 
numerically stable for all L. In contrast to Fig. [2j Fig. 
[3] demonstrates that tan 0^ stably converges with R by 
showing the numerically computed ratio of Wronskians 
in ( |16b[ ) for all values of R. These calculations use the 
same long-range potential as the calculations in Fig. [2] 
and include much larger values of L. 

The calculations of tan (f)i in Fig. [3] converge more 
rapidly as L increases. This trend is intuitive because, 
as L grows, the classical turning point at zero energy 
moves inward and X- becomes increasingly distinct from 
all other solutions. Moreover, these calculations of tan 0.; 
are stable out to very large R > 10^. This allows for 
a more consistent calculation of the MQDT parameters 
because the value of R at which 0i is actually deter- 
mined can equal the value of R that is necessary for a 
well-converged calculation of the MQDT parameters. Al- 
though either f or g could asymptotically coincide with 
X~ at zero energy, choosing g for this role defines our 
standardization because the MQDT parameters acquire 
very appealing qualities. The next section explores these 
qualities. 



zero-energy limit of Q, 



Q = - 



W{g,g) 
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is intimately related to the value of 
if (i 



7?— >oo 

(17) 

in (16b I. Clearly, 



in (17), the zero-energy limit of Q is the exact 
in (16b). Appendix [A| shows that 



same value as tan 



evaluating (16b) is stable independent of the phase 
given to the reference wave functions. Therefore, the 
zero-energy limit of Q is numerically stable. Fig. [4] shows 
the calculation of Q at 100 nK as a function of the zero- 
energy phase 4>i. 

For a given R^, Fig. |4] shows that there is a unique 
value of (pi (indicated by the green vertical solid line) for 
which Q vanishes at zero energy. This is the choice where 
g asymptotically coincides with x- at zero energy; there- 
fore, Q — Q &i zero energy defines our standardization. 
As evident from ^ and (10 1, the vanishing of Q guar- 
antees the maximal linear independence of / and g at 
long range. Hence, although a different standardization 
could guarantee that the zero-energy limit of Q is well 
behaved (e.g., the blue vertical dashed line in Fig. |4|, 
our choice of reference wave functions is the only choice 
that is maximally linearly independent at long range in 
the limit of zero energy. 

Applying our standardization. Fig. [5] illustrates the en- 
ergy dependence of the various MQDT parameters for the 
KRb model potential V^' = -Ce/R^ - Cs/R^ - Cio/R^°, 
where the dispersion coefficients are the same as in the 
previous section. The positive and negative energy pa- 
rameters are plotted together by defining the parameter 
JCi = \/\Ei\. The MQDT parameters are presented in the 
natural van der Waals units (3 = (2^C6//i^)^/^ = 143.9 
ao and Ep = 152.7 /zK. Each panel represents the result 
for a different partial wave L; note that a greater energy 
range is shown for higher L. In all cases the calculation 
is numerically stable, even in the threshold limit. These 
functions are smooth and hence easily interpolated. 

One striking feature, unique to our parameterization, is 
that all parameters vanish as powers of Ei in the Ei ^ 
limit. In this limit they are well approximated by simple 
analytic formulas. For alkali atoms, where Cs/R^ and 
Cio /R'^'^ make small corrections to C^/R^, these formulas 
can be derived using the —Cq/R^ potential alone. Their 
derivation is detailed in Appendix [B] and the results are 
summarized below. These parameters are conveniently 
parameterized in terms of a set of generalized, standard 
scattering lengths a^. 



IV. CALCULATING MQDT PARAMETERS 

Having specified 0^ in ([8|, the values of the MQDT 
parameters A, rj, Q, and 7 follow unambiguously from 



aL 



^2-(2L+3/2) 



2/(2L + l) 



(18) 



J(i/2 + 5/4)r(L + l/2), 
The threshold behavior of all four MQDT parameters is 
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IC, ip-') IC, (/3-') IC, {13-') 

(d)L = 3 (e)L = 4 (f)L = 5 



FIG. 5: (Color online) The MQDT parameters A (red solid curve), Q (blue dashed curve), rj (green dotted curve), and 7 
(black dashed-dotted curve) are shown for the KRb model potential V'" = -Ca/R^ - Cs/R^ - Cio/R'° with L = - 5. The 
parameter Jd = \/\Ei\ is defined merely for plotting purposes. 



given here, 



L+1/2 



' V ; V L y -r 32r(i + 7/2) ' 



19a) 

7 

(19b) 



fc2 



(2L + 3)(2L- 1) 



(19c) 




(2L+3)(2L-1) 



These formulas agree well with the numerical results for 
the exact long-range potential when Ei < 1. 



Expressions such as ( 18 1 have been derived before in 
the literature. For example, aq, which is the scatter- 
ing length of our reference wave function /, coincides 
with the semiclassical scattering length of Gribakin and 
Flambaum [57j for the —Cq/R^ potential. Likewise, us- 
ing an exact solution that was expressed using contin- 
ued fractions [12] [32], Gao performed a similar an- 
alytic treatment of the near-threshold MQDT parame- 
ters for the —Cq/R^ potential. To do so, he identified 
a set of standard constants CsL, Gao that are related to 
our equation (18 1 via asL,Gao = (o.l)^^'^^- While the 



of length. Moreover, our standard and universal refer- 



ence wave functions have a universal form (19b) for the 
corresponding phase shift. 

Gao conceives of a hierarchy of reference wave func- 
tions distinguished by a short-range quantum defect pa- 
rameter n'^ [13]. Our implementation of MQDT intro- 
duces alternative short-range phases 0i in (|8|. Our par- 
ticular choice for gives our reference wave function / 
a particular set of scattering lengths a^. Using equation 
(10) of |3D] and the zero-energy limit of equation (33) of 
[30j . this choice corresponds to. 



^= = 1/2 + L/4 



(20) 



Equation ( 20 ) gives an explicit connection between the 



analytic formulas of Gao and our formulation of MQDT 
for the pure —Ce/R^ potential. 

Moreover, a similar formulation of MQDT [T6| has de- 
rived the threshold behavior of three positive energy pa- 
rameters for a reference potential with arbitrary scatter- 
ing length a. These parameters are easily related to our 
MQDT parameters for the special case of L = 0. Ex- 
pressed in our notation, the results of jl6| read, 



77(L=o) > ~ak , 

fe^O 



treatments are equivalent, our parameters have units 



A(^L=o) - — > aok (1 + {a/ao - if) 
Q{L=o) -^^-^ 1 - a/ao • 



(21) 
(22) 
(23) 
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If a = ao here, these formulas are consistent with the 



threshold behavior of our set of MQDT parameters ( 19 1 



Although Appendix |B] explicitly derives the MQDT 
parameter threshold behavior for only the poten- 
tial, this analysis implies simple extensions of ( |19[ ) for 
any — potential. Here, n is any integer (n > 2). 
The zero-energy solution X- is well defined (and known 
analytically) for any potential of this kind. 



X- 



RJp 



{2L + l)/2v 



-s-cx R 



(24) 



where v — {2L + l)/(n — 2). Hence, for any potential 
asymptotically dominated by —Cn/R" and for all L, our 
standardization uniquely specifies the zero-energy limit 
of / and g, and the analysis of Appendix [b| is repeatable. 
As our standardization demands that all MQDT parame- 
ters go to zero in the limit i?i — 0, / and g have maximal 
linear independence at zero energy - even for this more 
general potential. A future report will consider the full 
theory for arbitrary n and L. 
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V. HIGH PARTIAL WAVE RESONANCES 

Performing a full, numerical calculation of scatter- 
ing observables allows the results of our formulation of 
MQDT to be tested. Here, the Johnson log-derivative 
propagator method [3T] produces numerically exact so- 
lutions to the coupled Schrodinger equations ([2]). This 
method propagates the log-derivative matrix Y to very 
long range R « 20, where it approaches a constant. After 
this full close-couphng (FCC) propagation, the asymp- 
totic limit of Y determines all scattering observables. 
The FCC calculation is known to be quite accurate, and 
it is our standard with which to compare the accuracy of 
MQDT. 

The prediction of high-L resonances in atomic colli- 
sions often requires many {N » 1) channels, causing 
the FCC calculations to become very time-consuming. 
The time required to compute the scattering matrix at 
a single energy and magnetic field is proportional to 
N^. Moreover, resonance widths decrease quickly with L. 
Thus, the prediction of resonance locations and widths 
requires many long FCC calculations. However, if the el- 
ements of K^'^ are weakly dependent on energy and mag- 
netic field, MQDT describes the same resonances after 
only a few iterations of a fraction of the full calculation. 
For instance, interpolating K^"^ on a coarse grid in mag- 
netic field greatly increases the numerical efficiency of 
calculating Fano-Feshbach resonances. 

Consider, for example, the collision of ^'^K -t- *^Rb 
in the lowest hyperfine states \Fk,Mpj^) \Fjn„Mpj^^) — 
|9/2, —9/2) |1, 1), where a number of Fano-Feshbach res- 
onances have been observed. Using the same model po- 
tential for both the FCC calculation and the calculation 
of K^'^ allows for a direct comparison between the FCC 
and MQDT methods. Our model adopts the accurate. 



-0.12' ' 
200 400 600 800 1000 

B{G) 

FIG. 6: The eigenphase shifts /xa are shown for the collision 
of ""K + *'^Rb over the range B = - 1000 G with a colhsion 
energy of 1 /iK. For this calculation the log-derivative matrix 
Y is matched to solutions in the long-range potential V^'' = 
-Ca/R^-Cs/R^-Cio/R^" at R^ = 45 ao, including channels 
with L = and L — 2. The two graphs merely display two 
different ranges of /i over the same range of B. 

short-range X and a ■^1]+ molecular potentials of 
[26] that were constructed by performing a global fit to 
the position of the L = Fano-Feshbach resonances. For 
consistency, our model also adopts the interaction pa- 
rameters of that describe the long-range forces. The 
van der Waals parameters Cq, Cg, and Cio describe the 
long-range dispersion forces; the electron-exchange inter- 
action is, 

Sex = Aexi?^"-e-'5-^ , (25) 

which is added to the triplet molecular state and sub- 
tracted from the singlet molecular state; and the dipole- 
dipole interaction is, 

(26) 

where a is the fine structure constant. 

Beyond i?m, the MQDT and FCC calculations involve 
slightly different Hamiltonians. On the one hand, the 
MQDT reference wave functions are solutions in the long- 
range potential of our choice. Because the anisotropic 
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dispersion forces dominate at long range, our choice is 

y"- = -Ce/i?^ - Cg/i?' - Cio/i?^" . (27) 

As a consequence, each channel differs only by a con- 
stant energy set by its hyperfine quantum numbers and 
subsequent Zeeman shift. This choice of V^^ ignores all 
other forces and all couplings between channels beyond 
i?ni- On the other hand, the FCC calculation considers 
the full Hamiltonian into the asymptotic region. 

However, Eex is vanishingly small beyond i? « 30 ao; 
hence, only H^^ is responsible for the difference between 
the MQDT and FCC calculations. H^d is very long range 
and creates a coupling between channels, but its inclu- 
sion beyond i?m = 45 ao makes a negligible contribution 
to the elastic cross section in this case. This allows for 
excellent agreement between the two calculations. In ap- 
plications where such weak longer-range couplings must 
be included, it should be straightforward to include them 
perturbatively, along the lines formulated, for instance, 
in [32] [35] [31]. 

If one chooses an where all channels are locally 
open, typically < 35 — 50 ao for alkali atoms, all 
resonant behavior is due to physics beyond this range 
and, therefore, approximately described by the MQDT 
parameters. Hence, the choice of Rra = 45 ao leads to a 
smooth that is easy to interpolate over a large range 
of collision energy and magnetic field. The eigenvalues of 
Ksr are known as eigenphase shifts ii\ = tanSx/n. Fig. 
[6] shows these eigenphase shifts as a function of magnetic 
field over the range B = - 1000 G. 

With this same choice of i?ni = 45 ao, a very coarse 
magnetic field grid of spacing 100 G allows for an accu- 
rate interpolation of K^'^, and MQDT accurately repro- 
duces the measured Fano-Feshbach resonances in L = 
— 2 states ^5] . This accuracy and the excellent agree- 
ment with the FCC calculation motivate using MQDT 
to quickly re-fit the singlet and triplet scattering lengths, 
producing our own accurate scattering model. Our fit 
includes the experimentally measured resonance posi- 
tions of [35] in all i < 2 states. It also includes the 
L = 2 resonance at 547.4(1) G reported in [35^ and 
recently confirmed by i37 . While retaining the value 
of Cq = 4.300 X 10^ atomic units, varying the scatter- 
ing lengths leads to a minimum reduced chi-squared be- 
tween the experimental and MQDT resonance positions. 
The optimal scattering lengths are as = —110.8 ao and 
at = -214.5 ao with x^cd = 0-83. 

Using this re-tuned Hamiltonian, our model predicts 
the position and width of the L = — 2 Fano-Feshbach 
resonances for "'"K -f- *^Rb collisions in their lowest hy- 
perfine states. For example. Fig. [7| shows the MQDT 
and FCC calculations for overlapping s- and d-wave res- 
onances. By only calculating Kg^. once every 100 G and 
interpolating over the range B = 0- 1000 G, MQDT re- 
produces the FCC calculation of resonance positions with 
an accuracy of < 1 mC. Moreover, since using MQDT to 
search for resonances only requires a magnetic field grid 
finer than the distance between any two resonances |38j . 




FIG. 7: (Color online) The elastic cross section for the col- 
lision of *°K -I- **''Rb with a colhsion energy of 1 is shown 
for the FCC calculation (black dots) , including channels with 
L = and L = 2. The FCC calculation is compared to the 
MQDT calculation (red curve) with Kst interpolated over a 
range of 1000 G. This curve is unchanged if the analytic for- 
mulas ( 19 1 are used instead of the numerical values for A, rj, 
and Q. 



the method allows enough numerical efficiency to ensure 
the discovery of all the Fano-Feshbach resonances in this 
range of magnetic field. 

Table [l] lists these resonance positions and widths. 
Even though some resonances are very narrow, finding 
the roots of det {Kqq + cot 7) determines quantitatively 
accurate resonance positions [10], where both K^'^ and 7 
are interpolated with ease. As our model calculates all 
experimentally measured resonances close to their posi- 
tions Bex, the unmeasured resonance positions of Table 
[l] are predictive with uncertainties on the order of cur- 
rent experimental uncertainties (< 1 G). Once our the- 
ory predicts a resonance at the position Bth, fitting the 
divergence of the scattering length near the resonance to 
the following form determines the width of an L = or 
L — 2 resonance [351. 



a{B) = flbg 1 



A 



B- B 



th 



(28) 



where Obg is the local background scattering length and A 
is the field width. Fitting the divergence of the scattering 
volume to the same form as equation ( 28 ) determines the 
width of an L = 1 resonance. 

For high-L resonances beyond L ~ 2, the resonance 
widths in K -|- Rb collisions become orders of magnitude 
more narrow. For example, MQDT predicts the widest of 
the L = 4 resonances to have a width < 1 /iG. The time 
required to perform a FCC calculation of resonances this 
narrow makes the comparison between the MQDT and 
FCC calculations challenging. However, predicting the 
position and width of high-i resonances remains simple 
within MQDT. Despite only performing a detailed, fully- 
coupled calculation on a magnetic field grid of spacing 
100 G, our method has found and characterized features 
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TABLE I: All of the Fano-Feshbach resonances in the range 
B = - 1000 G for the colhsion of ''"K + **^Rb in the state 
\Fk, MFk) \Fm, MFm) = |-9/2 - 9/2) |1, 1) are calculated 
using MQDT for a collision energy of 1 pK. These resonance 
positions Bth and field widths A are listed here with their 
associated partial wave quantum number L and compared 
with experimentally measured resonance positions Bex- Note: 
All magnetic field values are in units of Gauss. 
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18 orders of magnitude smaller than this. 



VI. CONCLUSION 

MQDT has been extended to high partial wave cold 
collisions. Our specific standardization of reference wave 
functions has produced a numerically stable calculation 
of high partial wave MQDT parameters that are smooth 
in energy and magnetic field. All of these parameters 
are described by simple power laws at ultracold ener- 
gies, and accurate expressions for these parameters in the 
threshold regime have been derived for potentials domi- 
nated by —Cq/R^ at long range. As an example, excel- 



lent agreement has been shown between MQDT and the 
FCC calculation of ultracold "^^K + ^''Rb scattering in 
their lowest hyperfine states. This calculation has also 
shown good agreement with experimental measurements 
of Fano-Feshbach resonances, and all of the L = — 2 
Fano-Feshbach resonances in the range of B = — 1000 
G have been reported. 
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Appendix A: The Instability of (j>i 

Here, we consider the calculation of the phase that 
defines the MQDT reference wave functions / and g. (f>i 
describes the particular linear combination of /(^.^o) ^iid 
5(0i=o) that coincides with a wave function of our choice 
at zero energy. We choose this wave function accord- 
ing to its asymptotic behavior and define our reference 
wave functions at i?x- Therefore, we can determine (pi by 
numerically propagating /(^.^o) ^-nd g(0.=o) from their 
boundary conditions at i?x to large R. However, numer- 
ical error causes this propagation to become unstable in 
the presence of a centrifugal barrier; hence, the numerical 
calculation of can also become unstable. 

For example, we consider the — 1/i?^ potential. In this 
potential, every zero-energy wave function is known in 
terms of the analytically known wave functions ( 15 ) at 



all R. These wave functions have well-known asymptotic 
behaviors and are exact solutions at small R. Therefore, 
their behavior at i?x determines 0.^ without any numerical 
propagation. Nevertheless, we can still attempt to deter- 
mine numerically by propagating the wave functions 
/(0.=o) and 3(0.^0) from their boundary conditions at R^ 
to large R. Then, studying the deviation of our numeri- 
cally determined value of (f)i from the analytically known 
value allows us to characterize the numerical instability of 
this calculation. Moreover, we identify a robust method 
for avoiding this instability that is easily generalized to 
any potential that falls off faster than 1/R^ asymptoti- 
cally. 

We rewrite the zero-energy wave functions f((t„=o) and 
g{(Pi=o) in terms of the analytic wave functions (15) by 
defining a constant 2x2 matrix C, 



/(0,=o) = ciix+ + ci2X- 
5(0,=o) = C21X+ + C22X- 



(Ala) 
(Alb) 



If we choose to let / coincide with x+ at zero energy, 
we can derive an expression for tan (pi in terms of the 



11 



elements of C by using equation (16a I, 



at large R, 



tan ( 



ciiW (x+, x+) + C12W X-) ^ C12 
C21W {x+,X+) + C22W {x+,X-) C22 



(A2) 

However, achieving this value of tan (pi numerically is 
not guaranteed. We can track numerical error by consid- 
ering the difference between a particular analytic solution 
and the same solution determined numerically. To this 
end, we expand the zero-energy wave functions x+ ^^id 
X- in powers of R at large i?. 



-R>>1 

x+ > 

R»l 

X- > 



L+l 



R 

R-^ 



BR 



L-3 



(A3a) 
(A3b) 



where i? is a known constant and the normalization of 
these wave functions is chosen such that W{x-,X+) = 
2L + 1. We then define two numerically determined wave 
functions x'+ and that have the same boundary condi- 
tions at i?x as the analytic solutions x+ and respec- 
tively. Moreover, we demand that x+ and x~ have the 
exact same normalization as their analytic counterparts 
asymptotically. 

Numerical error causes the wave functions x+ and x'_ 
to differ from the analytic solutions in two ways. First, 
only the leading order terms in their asymptotic expan- 
sions agree exactly. Hence, the asymptotic expansion of 
x'+ has a coefficient B' in its second highest order term 
that differs slightly from the coefficient B in equation 
Thus, we represent our error by the constant 
B — B' . Second, the numerical wave functions be- 
come slightly different linear combinations of the analytic 
solutions, such that x+ gains a contribution from x- that 
is proportional to the error S. differs from X- in anal- 
ogous ways, and we expand both of the numerical wave 
functions in powers of R at large i?, 



(A3a 



(A4a) 



X- 



-^K R-^ +6{R^+^ +BR^-^) . (A4b) 



The presence of numerical error also changes the wave 
functions f(^^=Q) and g(0.=o)- We call these numerically 
determined wave functions /^'^.^q-j and ff^^.^p)' they 
lead to the numerically determined phase We define 
/(0 =0) ^nd 5(0. by their boundary conditions at R^. 
Therefore, these functions are exactly the functions ( Al ) 



at i?x, but they take a slightly different form at large R. 
We approximate their large- i? behavior as the following. 



/(0,=O) > C11X+ + C12X- , 

R»l I , , 

5(,/,.=o) ^ C21X+ + C22X- ■ 



(A5) 
(A6) 



Using these wave functions in equation (16a) leads to an 
equation for tan0^ that depends on the product SR?^^^ 



tan^; ^ , (A7a) 

C2iM^(x+,x;)+C22M^(x+,x'-) 



.2L-3 



Ci2(2L + l) 



C2l4<5i?2i-3_C22(2L+1) 
fl-J-oo Cii 



> — for L > 1 and (5 7^ , 
C21 



(A7b) 
(A7c) 



For L > 1, the large- i? limit of tani/)^ approaches the 
wrong value C\\jc2\ if (5 is non-zero. Indeed, performing 
the actual numerical calculation produces this same value 
of tan^J. 

Since the value of tan <pi depends on the constants C12 
and C22 and these terms are dominated by numerical error 
at large i?, we deduce that finding the contribution to 
/(0i=o) and 5(0.^0) from x- at large R is numerically 
challenging when L > 1. In fact, letting either / or 5 
- at zero energy - coincide with any wave function with 
a contribution from x+ leads to an equation for tan (pi 
that depends on the constants C12 and C22, and the same 
numerical instability exists. However, if we instead let / 
or g coincide with x- at zero energy, using the numerical 
wave functions /(^.^g) and ff^^.^g) at large R leads to an 
equation for tan (j)'^ that reduces to the analytic value of 
tan (pi for all L. 

For example, if we let g coincide with x- at zero en- 
ergy, we can derive an expression for tan <pi by using the 
exact values of /(^^^g) and g(0.=o) in equation (16b), 



^ C2iVF(X-,X+) +C22M^(X-,X-) C21 

tan 0j = 7 —I- = . 

cnI^(x-,x;)+ci2VK(x-,x'-) cu 

(A8) 

Using the numerical wave functions Z^'^.^g^ and g^^.^g-j 
at large R leads to the following equations for tan (p'^ , 



, C21W {x-,x'+) + C22W {x-,x'-) . 
tan 0, = TTT-T r4 , (A9a) 



CiiW^(x-,x'+)+Ci2W^(x-,x'-) ' 



fl>>i^ C21 + C22S 

Cll + C12S 



C21 

Cll 



for all L and S << 1 . 



(A9b) 
(A9c) 



Here, tanc/)^ does not depend on R in the region i? >> 1, 
and tan</)^ approaches approximately the correct value 
even if the numerical error is finite. Of course an ac- 
curate value of tan (p'^ requires the numerical error to be 
small {S « C12 a nd S « C22), but the divergence seen 
in equation (A7b) does not appear. In this case, tan pi 



does not depend on the constants C12 and 022- Therefore, 
finding the contribution to /(0=o) and 5(0=o) from x- at 
large R is not necessary, and the numerical instability of 
calculating tan pi vanishes. 
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Appendix B: MQDT Threshold Behavior 

1. Introduction 

The MQDT parameters A, rj, Q, and 7 connect the 
reference wave functions / and g to well-known solutions 
in the limit i? — 00, as described in the equations ([o] 
- 111. By considering a simple long-range potential, we 



can represent the zero-energy limit of / and g in terms of 
analytically known zero-energy solutions. However, these 
solutions are inadequate to describe the large- i? behavior 
of / and g at non-zero energies, so we find a correction to 
the zero-energy wave functions via perturbation theory. 
With an accurate representation of / and g at small ener- 
gies in hand, we derive simple expressions for the MQDT 
parameter threshold behavior by matching these wave 
functions to either the energy-dependent wave functions 
/ and g or the function e^'^^. 



2. Zero-Energy Solutions 

For the simple long-range potential —Cq/R^, we can 
solve the Schrodinger equation, 



d^ 
cLR2 



L(L + 1)^ 

i?2 



i?6 



(Bl) 



analytically at £^ = 0. For all of Appendix [B| i? is in 
units of the natural length scale /3 — {2iJ,CQ/h'^y^^ of 
the potential ~Cq/R^, and E is in units of the natural 
energy scale Ep = /i^/2/Lt/3^, where /z is the re duce d mass. 
We describe particular solutions to equation (Bl ) at zero 



energy in terms of two linearly independent solutions x+ 
and X- defined by their asymptotic behavior, 

2L + l/2^L+l 

r(3/4-L/2) ' 

(B2a) 



X+ = 



= ^^-i(2L + l)(l/2i?' 



o-{L+l/2) D-L 

X- = VRJ.,,,^,,il/2R-) > r(X/2 + 5/4) ' 

(B2b) 

where J is the Bessel function of the first kind. 

For all energies, we define two linearly independent ref- 
erence wave functions / and g with the following bound- 
ary conditions at i?x << 1, 



f{R) 



1 



sin( / k{R')dR' + <j>) at i? = i?x , 
\/k{R) Jr^ 

(B3a) 

-1 

g{R) = —^= cos( / k{R')dR' + 0) aX R = R^ . 
^/k{R) Jr^ 

(B3b) 



Here, cf) is a phase that is constant in R a nd e nergv. and 
k = ^JE + l/i?6. The set of equations (B3) and their 

full radial derivatives define / and g. 



We then demand that g coincides (up to a normaliza- 
tion) with the solution x- zero energy. Hence, we 
rewrite / and g in terms of x+ and X- a-t zero energy by 
defining two constants of normalization A^i and A^2 and 
a constant phase a, 

j{E = ^)=N^ (X+ + tan ax- ) . (B4a) 
g(E = 0) = iViX- . (B4b) 

By considering the small- i? limit of our zero-energy solu- 
tions and reference wave functions. 



1 Ltt Stt 



R«\ 2 p3/2^- 

^+ ' '''' V 2i?2 4 

R<<1 2 „^/9 / 1 Lit Svr 



f{E ^ 0) R'/^ sin 

g{E = 0) 



^-i?3/2cos|- 



1 1 
2]^^2Rl^ 
1 1 



(B5a) 



(B5b) 



, (B5c) 



(B5d) 



we use the sets of equations ( B3 ) and ( B4 ) to determine 
the four unknown constants. 



2 ' 

1 Ltt Stt 

^ 1^ 8" 



tana = (-l)^+^sin 



2L + 1 



2sin(2iii^) ■ 



(B6a) 
(B6b) 

(B6c) 

(B6d) 



3. Perturbation Theory 

At zero energy, we know the wave functions / and g 
exactly; however, it is not immediately obvious whether 
or not the zero-energy wave functions are good approx- 
imations at large R >> 1 in the limit of i? — > 0. At 
small energies i? << 1, both / and g grow with R before 
reaching their asymptotic limits, but only / grows at ex- 
actly zero energy. Since / has a contribution from x+ 
and X- at zero energy, matching to finite-energy wave 
functions is straightforward. However, g is purely x- at 
zero energy, and the contribution to g from x+ at small 
energies is unknown. In order to match our zero-energy 
wave functions onto growing, finite-energy wave functions 
at large R, we must find the contribution to g from x+- 
This is accomplished by performing a perturbation in E. 

Because we plan to match wave functions at a finite i?, 
we choose a Green's function which preserves the bound- 



13 



ary conditions of g at i?x |39j . Hence, there is an integral equation for g at small ener- 

gies, 

Jo \iR<R' , 

I (jViATa)-! II n > n . 

(B7) 



.R 

g{R,E «1) ^ g{R,E ^0)+ G{R, R')g{R\ E = 0)AR! , (B8a) 

^0 

= N^x-{R) + ENiN2(^X+{R} x-{R')dR' -x-{R) x+{R')x-{R')dR'^ , (B8b) 

where we have used equation ( B4b ) to replace g at zero energy. We solve these integrals analytically and then expand 
them in powers of R at large R, 



R A 9-(2L + l) nl-2L 

(R')dR' ^ ^ - , — , , — , (B9a) 

' (2L + 3)(2L- l)7r (2L - l)r(5/4 + L/2)2 ' ^ ' 

' r , , , R»i 4cos Wivr) 2sin(^7r)i?2 

X-iR')x.iR>R' ^ (2^^3)(2L-1). + (2^:1). • (^^b) 

By using this correction to the zero-energy g and approximating / with its zero-energy limit, we have complete 
descriptions of the reference wave functions at large R in the limit of zero energy. 



f R»l /2i + l/2i?L+l 2-(^+l/2)i?-L\ 

/ > jVo — -, — -, + tan a — ; — -, 

■' E^o ^Vr(3/4-L/2) r(5/4 + L/2)y 

0-(L+l/2) p-L 

^ E^Q 'r(5/4 + L/2) 



(BlOa) 



EN^N2 



2L+5/2j^L+i 2-^+1/2 sin (^7r)i?-^+2 2-^+3/2 cos(^7r)i?-^ \ 

r(3/4 - L/2)(2i + 3)(2i - l)7r r(5/4 + L/2)(2i - l)7r r(5/4 + L/2)(2L -(- 3)(2i - l)7r ) ' 

(BlOb) 



4. Matching Wave Functions 

The MQDT parameters connect / and g with /, g, 
and e~'^^. Since we have analytic expressions for all of 
these wave functions at large R and small energies, we 
use equations ([o]- 111 to solve for the threshold behavior 
of the MQDT parameters. Moreover, by expanding these 
wave functions in powers of R and comparing like terms, 
we derive simple formulas. We can use the asympt otic ex- 



pansions of / and g at small energies in equations (BlOa 



and (BlOb), but we still need to find similar expansions 



of / and g in this same parameter regime; i? >> 1 and 
E^O. 

At large i?, we rewrite / and g in terms of spherical 
Bessel functions using Q and ([To]), 

/ UikR) cos 77 — n{kR) sin-q) 



kR 

— ^{n{kR) COST] + j{kR) sin T]) 



(Blla) 
(Bllb) 



Then, the small- argument expansions of the spherical 
Bessel functions unveil the behavior of / and g at large 
R and very small energies such that kR << 1, 



, 1 (fci?)^+i 

^i^TkW^-''''' 

1 (2L-1)!! 



1 (2i-l)!! 
^smr? 



Vk (kRY 



{kRf 

.^o^-y^"W^'^°''^l'+4i^2 

\L + 1 



1 {kRY 

7% (2i + 1; 



smry , 



(B12a) 



(B12b) 



where k — \/E. We see that / has a term proportional 
to and a term proportional to i?^^. Hence, we 

can compare this function with / term by term. The 
equation for / in ([o]), 

A-^'^f = J , (B13) 
yields two equations for the MQDT parameters A and ?7, 
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1/2 ^ r(3/4-£/2)fcWcos77 
Ar22^+i/2(2L + l)!! 



r(5/4 + i/2)(2L- 1)!! 



(B14a) 
(B14b) 



where matching powers of R^^^ leads to equation ( B14a 
and mat ching powers of leads to equation ( |B14b \ 



From (B14| we find that tan 77 oc Thus, for 



small fc, we use the small-angle approximation, sinry « 
J] and cosry 1, and define the generalized scattering 
length fli. 



^2-(2i+3/2) 



r(L/2 + 5/4)r(L + l/2) 



2/(2L+l) 



(B15) 



We then rewrite our expressions for A and rj in terms of 

= -(aifc)^+i/2 , (B16a) 
77 = (-l)^+i(aLfc)2^+i . (B16b) 

Here, we have used the relations (2n - 1)!! = 2"r(l/2 + 
n)/y^, where n is an integer, and sin (ttz) = 7''/r(l — 



z)r{z) with z = {2L + l)/4. For L > 1 and small k, we 
know the phase shift is dominated by a long-range phase 
shift proportional to k'^ 40j, but the derivation above 
only yields the short-range contribution because we are 
matching wave functions under the centrifugal barrier. 
Since we know the long-range contribution analytically, 
we simply correct our expression for rj by adding these 
contributions together. 



7? = (-l)^+^(aLfc)^^+^ 



37rr(£-3/2) ^^4 
32r(L + 7/2) 



(B17) 



The length scale helps to greatly reduce the amount 
of unnecessary constants in the derivation of the remain- 
ing MQDT parameters G and 7. Hence, we rewrite our 
wave functions in terms of ai,. At large i?, the zero- 
energy solutions become 



X- 



J?.>>i 1 





-{R/aL 



/a_L- 



(2L + 1)!! 
(2L- 1)!! 



(B18a) 
(B18b) 



and we have simple expressions for the zero-energy limits 
of / and g at large R, 



fl>>i. 



- jR/aL) 
{2L + 1) 

(2L- 1)!! 

(R/aLf 



L+l 



EH' 



{R/aL) 
{-l)^E 



E 



,L+1 



-{R/aL 
(2L+1)!! 



4L-2 (2L + 3)(2i-l)y (2i + 3)(2L-l 
We derive the threshold behavior of the MQDT parameter Q by using the equation for 5 in (|9]) 

A'/^g ^g + Gf . 



(B19a) 
(B19b) 

(B20) 



In this equation, we substitute g and g with their expansions in (B19b| and (B12b|, respectively, and replace A by its 
threshold value in equation (B16a). Thus, in the limit of large R and very small, positive energy such that kR << 1, 

(i?/ai)^ 
L)!! ' 

fc2 



we evaluate the LHS and RHS of equation (B20l separately, 



(2L- 1)!! 
^ ^""^ {R/aL)^ 

g + yj > \/aL 

'^^o {R/aLf 



1 



+ (-l)^(aifc) 



4L+2 



1 



4L - 2 

{kRf 

4L-2 ^ (2L + 3)(2L - 1) 



(2L+1)!! 



(B21a) 



(-l)^fc2 



+ {-irG]-v^ 



-{R/aL) 



L+l 



(2L-I-1)!! V(2i + 3)(2L- 1) 

(B21b) 



where E — k^. Hence, the first two terms on the LHS cancel exactly with the first two terms on the RHS, leading to 
the following equation, 



(-l)^(«Lfc) 



4L+2 (fi/«L)^+^ 

(2L + 1)!! 



(R/aL 

(2L + 1)!! V(2i + 3)(2L- 1) 



g )+(-!)'' L 



{R/Rl)"- V(2i + 3)(2L-1) 



(B22) 



In the limit of large R » 1, the _R^+^ terms on the RHS dominate the R ^ terms for all L and all fc, independent 
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of G'l therefore, neglecting the terms of order R ^ in this 
equation gives the threshold behavior of G, 



4L+2 



(2L- 



Jc2 

3)(2L 



1) 



(B23) 



We derive the threshold behavior of 7 in a way similar 
to the derivation of Q. Here, instead of matching to g, we 
need to match the small-energy limit of / and g to the 
function e""^ at large R. We again try to match wave 



functions at large R and very small energies such that 
kR << 1. Using equation (111, we define a constant of 
proportionality D, 



. ; , ^ -R>>1 „ — K i? 
tan7Ji + > De 



(B24) 



Then, using (B19) for / and g and using the kR << 1 
expansion of e^"", we arrive at the following equation. 



D 



{-kRY 



L+l 



n=0 



= -vaL 



-{R/aL) 
(2L + 1)!! 

(2L- 1)!! 

{R/aLf 



(2i + 3)(2L- 1) 
(nR)^ (-1)^k2 



tan 7 



4L-2 (2L + 3)(2L-1) 



+ (-l)^tan7 



(B25) 



where E = —k^. 

For i = 0, we take the expansion of e^'^^ out to first 
order in k (n = 1) and neglect terms of order k^. 



D{1 — kR) = \/ao (1 + tan7 — tan7i?/ao) 



(B26) 



Matching constant terms and terms of order R leads to 
the following two equations with two unknowns, 



D 

-Dk 



Hence, 



D 

tan 7 



/oj^ (1 + tan 7) 
tan 7 



\/aoK 
1 

1 — aoK 



(B27a) 
(B27b) 



(B27c) 
(B27d) 



and we have a simple formula for tan 7 in the threshold 
limit K << 1 when L = 0, 



tan 7 ~ aoK for L — 



(B28) 



For L > 0, we immediately see that matching powers 
of R in equation (B25 ) is problematic due to the terms of 
order on the RHS. Therefore, instead of matching in 
the limit kR << 1, we simply let the wave functions take 
their asymptotic forms as i? — 00, where kR >> 1 even 
though K << 1. That is, e~'^^ — > 0, and / and g are still 
well approximated by their E —i' limits in (B19|. As R 



becomes very large, these wave functions are dominated 
by their contributions from Even the term in g 

proportional to R^^~^^ is dominated by i?^"*"^ for L > 0, 
and equation (B24) takes a simple form. 



L + l 



= 



-(R/ol) 

(2L + 1)!! V (2L + 3)(2L-1) 



tan 7 



(B29) 

Hence, we have the following simple formulas for the 
threshold behavior of 7, 



aoK 

^ (2L+3)(2L-1) 



for L = 
for L > 



(B30) 
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